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retical uncertainty in the determination of the critical value of Mh 
for vacuum stability to 1 GeV. While A at the Planck scale is re- 
markably close to zero, absolute stability of the Higgs potential is 
excluded at 98% C.L. for Mh < 126 GeV. Possible consequences of 
the near vanishing of A at the Planck scale, including speculations 
about the role of the Higgs field during inflation, are discussed. 
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Abstract 



1 Introduction 



The value of the Higgs mass (Mh) favored by present ATLAS and CMS data [1,2], Mh ~ 125- 
126 GeV, is intriguing: it is quite close to the minimum Mh value that ensures absolute 
vacuum stability within the Standard Model (SM) which, in turn, implies a vanishing Higgs 
quartic coupling (A) around the Planck scale. In order to assess if the measured Higgs mass 
is compatible with such a peculiar condition, a precise computation is needed. 

The state-of-the-art analyses of SM vacuum stability before the latest LHC data were 
done at the next-to-leading order (NLO) level [3-10]. This is based on two-loop renorma- 
lization-group (RG) equations, one-loop threshold corrections at the electroweak scale (pos- 
sibly improved with two-loop terms in the case of pure QCD corrections), and one- loop 
improved effective potential (see [11] for a numerically updated analysis). 

With this paper all the ingredients necessary for a complete next-to-next-to-leading order 
(NNLO) analysis in the strong, top Yukawa and Higgs quartic couplings become available. 
In particular, complete three-loop beta functions for all the SM gauge couplings have been 
presented in [12], while the leading three-loop terms in the RG evolution of A, the top Yukawa 
coupling (y t ) and the Higgs anomalous dimension have been computed in [13]. However, as 
pointed out in [11], the most important missing NNLO piece for the vacuum stability analysis 
are the two-loop threshold corrections to A at the weak scale due to QCD and top Yukawa 
interactions, because such couplings are sizable at low energy. The calculation of such terms 
is presented in this work. 

The relation that connects A to the Higgs mass and to the the Fermi coupling (C? M ) can 
be written as 

C M 2 

A(a*) = -^+AA(a*), (1) 

where AA(/i) denotes the sizable threshold corrections arising beyond the tree level. Given 
the rapid variation of A around the weak scale (see fig. 1), these corrections play a significant 
role in determining the evolution of A up to high energies. Computing AA(/i) at the one 
loop level, using two-loop beta functions for all the SM couplings, and varying the low- 
energy matching scale between M t /2 and 2M t , leads to a ±2 GeV error on Mh [11]. The 
NNLO finite terms that we compute here allow us to reduce this error down to ±0.7 GeV. 
While this work was in progress an independent calculation of the two-loop Yukawa-QCD 
contributions to AA(/i) has appeared [14]. Our result agrees with the one in ref. [14] for 
these contributions. However, our analysis includes also the two-loop terms coming from 
the Yukawa sector and can be considered the first complete NNLO evaluation of AA(/i). 
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Figure 1: Left: SM RG evolution of the gauge couplings g\ = y / 5/3g', <?2 = 9, 53 = 9s, of the 
top and bottom Yukawa couplings (yt,Ub)> an d of the Higgs quartic coupling A. All couplings are 
defined in the MS scheme. The thickness indicates the drier uncertainty. Right: RG evolution of 
A varying Mt and a s by ±3o~. 

We stress that both these two-loop terms are needed to match the sizable two-loop scale 
dependence of A around the weak scale, caused by the —?>2yfg 2 s + 30yf terms in its beta 
function. As a result of this improved determination of AA(/x), we are able to obtain a 
significant reduction of the theoretical error on compared to previous works. 

Putting all the NNLO ingredients together, we estimate an overall theory error on Mh of 
±1.0 GeV (see section 3). Our final results for the condition of absolute stability up to the 
Planck scale is 

r ^ , (M t [GeVl - 173.1 \ (oi s (M z ) - 0.1184 \ 
M h [GeV] > 129.4 + 1.4 ( t[ ^ J - 0.5 I |^ J ± 1.0 th . (2) 

Combining in quadrature the theoretical uncertainty with the experimental errors on M t and 
a s we get 

M h > 129.4 ± 1.8 GeV. (3) 

From this result we conclude that vacuum stability of the SM up to the Planck scale is 
excluded at 2a (98% C.L. one sided) for M h < 126 GeV. 

Although the central values of Higgs and top masses do not favor a scenario with a 
vanishing Higgs self coupling at the Planck scale (M P i) — a possibility originally proposed 
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in ref. [15] and discussed more recently in ref. [5,6,8, 16, 17] — the smallness of A around Mpi 
is quite remarkable (see fig. 1). Motivated by this observation, we have explored in more 
detail the robustness of the predictions for and M t assuming special boundary conditions 
on A and its beta function around Mpi, as advocated in [16]. We have also critically examined 
scenarios where the Higgs field plays a role during inflation. This could happen because of 
a non-minimal coupling to gravity that flattens the SM potential close to Mpi [8], provided 
A(M P1 ) is positive. Alternatively, the Higgs field could have caused inflation while it was 
trapped in a second unstable minimum of the potential that appears near Mpi if A(Mpi) is 
positive and very close to zero [6], provided some non-SM mechanism is introduced to exit 
inflation [18, 19]. Beside some technical problems, these frameworks could become viable 
only if Mh satisfies the stability condition in eq. (2). We therefore conclude that both 
these possibilities are not favored by present data, unless M t is below about 172 GeV or 
new-physics threshold corrections at the high scale modify the shape of the SM potential. 

The paper is organized as follows: in section 2 we present the calculation of the threshold 
correction AA(yu). The numerical results for the condition of vacuum stability and, more 
generally, for the structure of the SM Higgs potential up to very high field values are discussed 
in section 3. The implications of these results for Planck scale physics are discussed in 
section 4. The results are summarized in the conclusions. We include also an appendix in 
which a ready-to-use expression for the two-loop effective potential is presented. 

2 The two-loop threshold correction to A(/i) 

In this section we present our main new result, namely the calculation of the two-loop 
contribution to AA(/i). We first obtain the y\ and y\g 2 s terms from the calculation of the 
Higgs mass via the effective potential. Then, we present the full result for the two-loop QCD 
and Yukawa contribution to AA(/i) in the SM with the electroweak gauge couplings switched 
off (the so-called gauge-less limit). 

2.1 Two-loop corrected Higgs mass from the effective potential 

We write the SM potential for the Higgs doublet H in the usual way: 
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so that, up to negligible width effects, the pole Higgs mass Mh is the solution of the pole 
equation 

M 2 h = -m 2 + Z\v 2 + U hh (M 2 h ) , (5) 

where m 2 , A and v are MS renormalized quantities and Uhh{p 2 ) is the Higgs self-energy 
(two-point) function, with external four- momentum p. We rewrite this equation as 



Ml 



[_ m 2 + 3X v 2 + n fch (o)] + [n hh (M 2 h ) - u hh (o)] 

[M 2 h ] v + ATl hh (Ml) 



(6) 



This step is convenient because the last term (which is computationally challenging) only 
gives corrections suppressed by the small Higgs coupling, in view of the smallness of Ml = 
2Xv 2 at tree level. The first dominant term can be expressed in term of derivatives of the 
effective potential, V e s. Writing the effective potential as a sum of the tree- level part V plus 
radiative corrections AV 

V eS = ~h 2 + h 4 + AV , (7) 

one finds 

[Ml\ 



d 2 V P 



cff 



I V 



(dhy 



h=v 



where v is the h vev at the minimum of the effective potential, determined by the minimiza- 
tion condition 



dV e 



cff 



dh 



h=v 



-m 2 h + A/T 



dAV 
dh 



(9) 



- h=v 



As usual, it is convenient to consider m 2 as a free parameter fixed in terms of v by the above 
equation, arriving at 

1 dAV d 2 AV~ 



Ml 

Defining the operator V"^ as 1 



2Xv 2 



h dh 



+ 



(dhy 



(10) 



h=v 



V 



1 d 



+ 



d 2 



hdh (dh) 



h=v 



and noting that 2\v 2 = V^Vq, we can simply write 
form for the Higgs mass: 

V 2 m V eS + AU hh (Ml) 



111 



V^Vee, obtaining the following 



Ml 



(12) 



It gives the Higgs mass squared as the sum of two terms. The first is the Higgs mass 
obtained from the potential; this is not the complete pole Higgs mass and must be corrected 



1 Notice that the term in linear in field-derivatives automatically takes into account the cancellation 
of /i-tadpoles (or alternatively, the minimization condition to get the right v). 
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for nonzero external momentum effects, which are taken care of by the last term, AHhh(M^). 
It is a straightforward exercise to verify that this expression for pole mass is independent of 
the renormalization scale fi. In particular, one can easily prove that 

d^Mv = -HMl] v , (13) 

^-AH^M 2 ) = 2 7 [Ml - AU hh (M 2 )] , (14) 

where 7 is the Higgs anomalous dimension, describing its wave-function renormalization, 
7 = dlnh/dhi //. 

Using eq. (12) and the one-loop result for V e s in eq. (74) of the appendix one obtains the 
one-loop Higgs mass correction. The explicit one-loop result for the pole mass is 

Ml = 2Xv 2 + , (15) 

with 

5iM 2 h = {3y 2 (4m 2 - M 2 h )B (m t , m u M h ) + 6AV(34 - 6 + ttv^) 

v 2 

--(3<7 4 - 8A<? 2 + 16X 2 )B (m w , m w , M h ) 
-^(3G 4 - 8AG 2 + 16X 2 )B (m z ,m z ,M h ) 

o 

+2m 2 w [g 2 - 2X(£ W - 1)] + m| [G 2 - 2X{t z - 1)] } , (16) 

where G 2 = g 2 + g' 2 . All parameters on the right-hand side (including v) are MS running 
parameters (with the exception of M 2 , which appears through the external momentum 
dependence of the Higgs self-energy). As eq. (74) was obtained in the Landau gauge, v in 
eq. (16) represents the gauge and scale-dependent vacuum expectation value of the Higgs field 
as computed in the Landau gauge. Similarly the AH/ l / l (M|) contribution in that equation is 
computed in the Landau gauge. In eq. (16) 

. r 1 , (1 - x)m 2 + xm? - x(l - x)m 2 - ie , 
B (m a ,m b ,m c ) = - / In ^ b — -± '—^ dx , 17 

Jo V 
and £ x = ln(m 2 //z 2 ), with m x the running mass for particle x (m t = y t v/\/2). One can 
explicitly check, using the RGEs for these parameters, that this expression for M 2 is indeed 
scale-independent at one-loop order. 

Neglecting gauge couplings and setting M% = 2Xv 2 in the one-loop terms, one obtains 
the approximate expression 

5 1 M 2 h ~ ^ [A(2 + 3£ t ) - 3y 2 e t ] . (18) 
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Figure 2: Two-loop vacuum diagrams that give the dominant contribution (depending only on g s 
and yt) to the SM two-loop effective potential. 



To compute eq. (12) at the two-loop level one can use the two-loop effective potential [20] 
to calculate [M^] y and the general results for two-loop scalar self-energies in [21] (supple- 
mented by the results on two-loop momentum integrals of [22]) to calculate AII/^M^). If 
we only keep the leading two- loop corrections to M\ proportional to yf, yfg 2 , dropping all 
subleading terms that depend on the electroweak gauge couplings or A, our task is simplified 
dramatically. First, in the two- loop effective potential we only have to consider the diagrams 
depicted in fig. 2. Their contribution can be extracted from the expressions for Vy and Vpv 
in the appendix. Second, in the two-loop term An^(M^) we can substitute the tree-level 
2Xv 2 , so that 



value M\ 



AnS(M, 2 ) 



n$(2A« a )-n$(o). 



(2), 



(19) 



It is then clear that the two-loop contributions coming from that term are proportional to 
A and are therefore subdominant. In this section we nee: lect AngJ(M^) completely. 



To find the expression for the Higgs mass at two-loop precision, we must also take into 
account that M\ has to be evaluated with one-loop precision in the argument of the one-loop 
term An^(M^). Putting together all these pieces, keeping only the y\ and y\g 2 s terms, we 
arrive at the following two-loop correction to eq. (15): 



5 2 M 2 h 



(4tt) 4 



7T 



2\ 1 



(20) 



The expression for as a function of A can be inverted to obtain A(/i) as a function 
of the pole Higgs mass M h . To express A(/z) in terms of physical quantities (i.e. and 
the pole masses Mz, M\y, and M t ) the relations between physical and MS parameters are 
needed. At the level of accuracy we are working only the relation between the yt((i) and M t 
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and the one between v(n) and G M are required. They are given by: 



y 2 M = 2V2G.M? 
1 



v\n) 



+ 



V2G, (4tt)2 [ 



1 + lj4^ 9 ° {3LT " 4) + Wf^ G » M * { ~ 9Lt + n] 



3M 2 (2L T - 1) + M^(5 - 6L W ) + -Mf (5 - 6L Z ) 



where Lx = ln(M^//i ), with masses in capital letters denoting pole masses. 



(21) 



(22) 



We find: 



A(/i) = ^ + A«M+A< 2 >(/i) 



(23) 



with 



A (2) (A*) = Ag, D , Iead . (/i) + Ag kiIead . (/.) + ..., (24) 

where the ellipsis stands for the subleading terms neglected in this section. The known one 
loop term is 

6(L H - L w )Ml 



A«M 



y 2 i 

2 ^vr) 2 



M 2 -M^ 



- 8 (2M^ + M|) - 2(-3 + QL T )M 2 M 2 



+M 4 (l9 - 15L H + 6L W - 3V3nj + 12(M 2 - AM 2 )M 2 B {M U M t , M h ) 
+2 (M 4 - AM 2 h M^ + 12M^) 5 (M W , M w , M h ) 
+ (MjJ - 4M 2 M| + 12M*) £ (M Z , M z , M & ) 

+M 2 



6M 2 M 2 



2(8L W - 7)M^ + (8L Z - 7)M| - "T z ^ 



Mf-M 2 , 



(25) 



and the leading two loop QCD and Yukawa terms are 



A 



(2) 

QCD.lead 



A 



(2) 

Yukjead 



(A*) 



G 2 M 4 



(4tt) 



64^ 2 (/i) (-4 - 6L T + 3L 2 .) 



(4vr) 4 



(30 + 7T 2 + 36L r - 45L 2 .) 



(26) 
(27) 



The above expression for A(/i) has the correct dependence on the renormalization scale /z, 
so that both sides of (23) evolve with /i in the same way to the order we work. 



2.2 Two-loop contribution to AX(fi) in the gauge-less SM 

In this section we go beyond the leading yf and yfg 2 contributions to AA(/i) discussed 
in the previous section and derive the QCD and Yukawa corrections in the SM with the 
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electroweak gauge couplings switched off. We first discuss the two-loop renormalization of 
the Higgs sector of the SM and then derive the 2-loop relation between A(/x) and the physical 
parameters G^,M t , and M^. In our derivation we follow closely ref. [23]. 

We start from the unrenormalized Higgs potential, eq. (4) written in terms of bare quan- 
tities, and set m 2 — > m 2 — 5m 2 , X — > X r — 5\, v — > v r — Sv. Then, assuming 5m 2 , 5X, 5v of 
0(a), we obtain 



V = V r - 5V , 



(2f 



where, putting m 2 = X r v 2 , 



Vr 



X, 



2\2 



G + G- (G + G- + h 2 + G ) + - (h 2 + G 2 ) 



+ X r v r h[h 2 + G 2 + 2 G + G-] + -M 2 h h 2 , 



(29) 



with Mh = 2X r v 2 , while, up to two-loop terms, 



5V 



5X 

+ 



2\2 



G + G- (G + G- + h 2 + G ) + - (h 2 + G 2 ) 



X r 



Sv 2 
2v, 



+ 



(5v 



2\2 



+ v r 5X[ 1 - 



1 



2^2 



SM 2 h 2 + v r 8r 1 



h [h 2 + G 2 + 2G + G-] 
5v 2 ' 



In eq. (30) 



SM 2 h 



8i 



X r 5v 2 + v 2 5X 1 



5v' 



2v 2 



5m 



h 



X r 5v 2 + vl5\ 1 



5v" 



5m z 



(30) 

(31) 
(32) 



and Sv 2 is related to Sv through Wv 2 — 5v 2 = v r — 5v. 

Following ref. [23] we require the cancellation of the tadpole contribution by setting 

5v 2 \ T 



5t 1 



2vi 



Vr 



(33) 



where iT is the sum of the tadpole diagrams with the external leg extracted. We identify 
M 2 in V r with the on-shell Higgs mass leading to the condition 



8M 2 = ReU hh {M, 



h) i 



(34) 



where the H/ l / l (M|) in the above equation includes only the contribution of the self-energy 
diagrams because eq. (33) is enforced. A third condition can be obtained by looking directly 
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at the muon-decay process. At the two-loop level we can write 



V2 H 



A WW 2 Mw 

IV + M Wq B w + 



A ww Vw 



(35) 



where t> is the unrenormalized vacuum, A ww = A WW {Q) is the W self-energy evaluated at 
zero external momenta, Vw and Bw are the relevant vertex and box contributions in the fi- 
decay process and Mw is the unrenormalized W mass. Performing the shifts v q — > v 2 — 5v 2 , 
M Wo — > M w — SM W , where SM W = ReA W w(M w ), and working at the two-loop level we 
arrive at 



1 



1 



.4 



WW „ , A WW 5M W 



V2G U V2G U 1 M w 



E + 



M w 



Aww 



A WW V W i r /\ /t2 r> I c 2 

+ — ^72 — + 5M w B w >+Sv 
M w 



(36) 

where E = Vw + M w Bw- We identify the renormalized vacuum by v 2 = l/(v / 2G M ), then 
5v 2 is defined to cancel the contribution of the curly bracket in eq. (36). 



Our choice of renormalization conditions implies that the renormalized quartic Higgs 
coupling is set equal to 

(37) 



while eqs. (31-36) can be used to obtain the correction 5X. Writing 
where the superscript indicates the loop order, we have 



(3* 



5A« 
6\W 



+ 



+ 



A 





{ A {1) 
^ww 

{ M w 






rt-WW 

{ M w 




^■ww 
M 2 W 




( 4 (1) 

^WW 


v ww° 1V1 W 


( A W 

^WW 


M w \ 


\ M w 



M 2 h 



M 2 h 



Ren«(M, 2 ) + 



ReIT (2) '- 2 



- - 



L hh 



Ml 



(Ml) + 



v r 

V r 



(39) 



Ren«(M 2 ) + 



3 T« 



2 V r 



|W T/(l) 



^ w Vy , r(l) M 2 
Jl V 



(40) 



The connection between A r , as defined in eq. (37), and X(/i) can be easily derived using 

A r -8X= A(/i) - 5\ , (41) 
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or 



(42) 



In eqs. (41)-(42) 5\ is the counterterm associated to A(/x), i.e. the counterterm that subtracts 
only the terms proportional to powers of 1/e and 7 — ln(47r) in dimensional regularization, 
with d = 4 — 2e being the dimension of space-time. Concerning the structure of the 1/e 
poles in 5X and 8X, one notices that it should be identical once the poles in 5X are expressed 
in terms of MS quantities. Then, after this operation is performed, a finite A(/x) is obtained. 



Specializing the above discussion to the two-loop case we have 

a 



am 



V2 



M 2 h - 8\V\^-8\ { -% n + A 



from which we identify the one- and two-loop contributions entering eq. (23), 

A«(/i) = -*A< l V, 
which reproduces the one-loop result of ref. [23] and 



-5A {2) L+A 



(43) 



(44) 



(45) 



In eqs. (43)-(45) the subscript 'fin' denotes the finite part of the quantity involved and A is 
the two-loop finite contribution that is obtained when the OS parameters entering the 1/e 
pole in 8a^' are expressed in terms of MS quantities, the finite contribution coming from 
the 0(e) part of the shifts. 

Differentiating eq. (43) with respect to fi, the known two-loop beta function for the Higgs 
quartic coupling is recovered. It should be recalled that the right-hand side of eq. (43) is 
expressed in terms of physical quantities, then the dependence on /i in that equation is 
explicit. To obtain the correct two-loop beta function, one has first to differentiate with 
respect to // and then to express the one-loop part in terms of MS quantities. 

The computation of A^(/i) in the full SM is quite cumbersome, see eq. (40). However, 
the calculation can be greatly simplified if one considers the gauge-less limit of the SM in 
which the electroweak interactions are neglected, i.e. the gauge couplings g and g' are set 
equal to zero. In this limit, eq. (40) simplifies to 



(5\W - A) 



g.l. 



+ 




RenS(M^) + 



T(2) 



RenS(M 2 



M 2 W 



2 v r 



A 



g.l. 



(46) 



g.l. 
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where the subscript g.l. means that we have considered in the various self-energies only 
diagrams involving the top and bottom quarks, the Higgs and the Goldstone bosons, the 
latter with vanishing mass, and the limit g,g' — >■ is taken. 

Using eq. (46) we compute the QCD and the Yukawa contribution to A*- 2 ^). The top 

(2) 

Yukawa-QCD contribution, Aqq D (/x), is obtained evaluating the relevant diagrams via a 
Taylor series in Xht = M 2 /M 2 up to fourth order 



V QCD 



GlMf 



(4vr) z 
+ x ht ( 35 



16 (-4 - 6L T + 3L 2 T ) 



2n 2 



12L T -12L 2 T )+x 2 ht — + x ht 



61 
135 



3 1223 4 43123 

+ x ht - 



3 1 V ~" r 135 ' ~" r 6300 ' ™1323000 

where N c and Cp are color factors (N c = 3, Cp = 4/3). Equation (47) shows that the 
series converges very fast. Our result is in agreement with ref. [14], the numerical difference 
between eq. (47) and the expression of ref. [14] for Mh ~ 125 GeV being negligible. 

The Yukawa contribution, A Yi ; k (/i), is (neglecting the small bottom Yukawa) 

V2G%Mf 



^Yuk^) 



n: 



lQB (M t ,M t ,M h )(-l + 2L T ) 



+N C 



+ x ht ((1 + AB (M t , M t , M h ) - 2L T ){\ - 2L T )) 

16 + j^vr 2 + 32B (M t , M h , M t )(l + 2L T ) - A8L 
(929 16 

-x ht [—+ — 7T 2 + A8B (M h , M h , M h ) - If 



^-L T + 58L 2 T 



129 16 

— + — vr 2 + 48B (M h , M h , M h ) - 16L H (1 - l t ) 
\ 6 3 

+B (M t , M h , M t ) ( y + 32L T ^j + ^ 
+ 4t + \^ - \ L h + B (M h , M h , M h ) (27 - 18L 

+10L T L H + 12L 2 T + B (M t , M h , M t ) + 4L T ^j ^ 



■ T ) + A0L T 



» /1181 7T 2 61 , Ti jt \ 59 r 

+ xt 1 B (M h , M h , M h ) H L 

'" 1 900 2 30 0V ' ' ; 90 



3 

90 £ 



-A So (M t ,M ft ,M t )-§ 
35 63 



131 2 /729 135 . 

vr 2 + V3 

6 V 2 ! 

/-225V3 

+ 71 I 



j-JT 

63 

»-^V3n)s 2 -lllL H + 36L 2 H 
4 / 

^ + 18V3L^ + 75 + 4 72C3 ]} , 
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(48 



where B (M h ,M h ,M h ) = 2 - L H - n/V3 and S 2 = 4/(9 y/S) Cl 2 (7r/3) = 0.260434138... 
In eq. (48) the terms proportional to AT? and AT? were computed exactly while the ones 
proportional to N c were computed via an asymptotic expansion in the large top mass up to 
x\ t terms exploiting the asymptotic-expansion techniques developed in ref. [24]. The part 
independent of N c in eq. (48) was computed using the results for the two-loop on-shell master 
integrals of ref. [25]. 

We end this section by commenting on the size of the terms suppressed by powers of 
Xht with respect to the yfg^ and yf contribution in eqs. (47, 48). While in the QCD case, 
eq. (47), the Xht suppressed terms are indeed smaller than the yfg^ contribution, the same 
is not true in the Yukawa case, eq. (48), where the x^t terms are actually larger than the y\ 
contribution. 



3 Extrapolating the SM up to the Planck scale 

A full NNLO computation of the Higgs potential requires three main ingredients: 1) the two- 
loop effective potential; 2) three-loop beta functions for all the relevant couplings; 3) two- 
loop matching conditions to determine the initial values of the couplings at the electroweak 
scale. As anticipated in the introduction, all these ingredient are now available for the QCD, 
Yukawa and Higgs quartic couplings. In this section we first discuss the structure of the 
two-loop potential and the numerical inputs at the electroweak scale, and then present the 
final numerical results for the stability condition in the M^-Mt plane. 

3.1 The two-loop effective potential 

The SM effective potential is known up to two-loops [20] . Its explicit form in a ready-to-use 
expression is given in the appendix. For large field values (h v), the potential is very well 
approximated by its RG-improved tree-level expression, 

VT(h) = , (49) 

with fi = 0(h). For this reason, if we are interested only in the condition of absolute 
stability of the potential, we could simply study the RG evolution of A imposing the condition 
A (A) > for any value A up to the Planck scale (as for instance done in [4]). Given 
that A reaches its minimum value before Mpi, independently of its initial condition at the 
electroweak scale, the minimum Higgs mass ensuring vacuum stability corresponds to the 
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initial value of A such that at some scale Ao 

A(A ) = /3 A (Ao) = , 



d 



(iln/i 



AO*) 



(50) 



This is indeed the condition analyzed in ref. [14]. In principle, a more accurate determination 
of the minimal Mh ensuring vacuum stability is obtained taking into account the full structure 
of the Higgs potential at the two-loop level. In practice, the determination of Mh obtained 
by the condition (50) differs by about 0.1 GeV from the one determined by the absolute 
stability of the RG-improved two-loop potential. 

In the following we are interested also in analyzing the shape of the Higgs potential close 
to the Planck scale and in the scale where the instability occurs (as a function of and 
M t ). To this purpose, the study of the RG evolution of A is not sufficient and the complete 
structure of the effective potential at the two-loop level plays a significant role. As pointed 
out in [3], one can always define an effective coupling X e s(h) such that for h 3> v the two- loop 
effective potential assumes the form 

Aeff(^) 



V eS (h) 



-h 4 



(51) 



The explicit two-loop result for X e g(h) can be easily obtained from the two-loop potential 
and is given in the appendix. We report here the simplified expression obtained when, in 
the two-loop term, we take into account only the contributions from the strong and the top 
Yukawa couplings 2 [3]: 



Xeff(h) 



X(h) 



(4tt) 



(52) 



+ 



(47T) 



8g 2 s (3rf - 8r t + 9) - -yf ^ 3r 2 - 16r t + 23 , 



7T 



Here all couplings are evaluated at the scale determined by the field value (/i = h), the index 
p runs over particle species, N p counts degrees of freedom (with a minus sign for fermions), 
the field-dependent mass squared of species p is rrvi{h) — + n p h 2 and C p is a constant. 
The values of {N p , C p , fi p , k p } within the SM are: 



(53) 



2 At high scales, the electroweak gauge couplings g' and g become comparable in size to y t and g s (see 
fig. 1), but their contribution to \ e s(h) turns out to be numerically small so that eq. (52) is a very good 
approximation. 
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Figure 3: Evolution of the Higgs coupling A(/i) and its beta function, eq. (50), as a function of the 
renormalization scale, compared to the evolution of the effective coupling A c g(/i), defined in eq. (51), 
as a function of the field value. Left: curves plotted for the best-fit value of Mt- Right: curves 
plotted for the lower value of Mt that corresponds to \{M-p\) = 0. 

The factor 

V(h) = I 70)dln/i , (54) 

where 7 = din h/ din fi is the Higgs field anomalous dimension, takes into account the wave- 
function renormalization. We have also defined r p = ln[/t p e 2r( -' , '- ) ]. 

The difference X e s(h) — X(h) is positive, as illustrated in fig. 3. As a result [3], at a 
given field value the potential is more stable than what guessed from the naive expectation 
based on the RG-improved tree- level potential in eq. (49), with fi = h. We finally notice 
that the difference \ e g(h) — X(h) gets suppressed at large field values, especially when A 
reaches its minimum close to the Planck scale. This is expected according to the following 
two observations: 1) the difference between A e g and A can be reabsorbed by a shift in the 
scales at which the two couplings are evaluated, up to finite two-loop corrections; 2) this 
shift has a small impact at large field values given the corresponding vanishing of (3\ (see 
fig. 3). 
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3.2 Inputs at the electroweak scale and threshold corrections 



As far as the SM gauge couplings are concerned, we can directly use results in the literature 
for the couplings in the MS scheme. In particular, from a global fit of electroweak precision 
data, performed with the additional input M h ps 125 GeV, the following MS values of the 
electromagnetic coupling and the weak angle renormalized at Mz are obtained [26]: 

ai£ = 127.937 ± 0.015 , sin 2 6 W = 0.23126 ± 0.00005 . (55) 

From these we derive 

a 2 \M z ) = a-^sin 2 #w = 29.587 ±0.008 , (56) 

a Y 1 (M z ) = cos 2 6 W = 98.35 ± 0.013 . (57) 
For the strong coupling we adopt 

a s (M z ) = 0.1184 ±0.0007 [27] (58) 

such that, including 3 loop RG running up to M t and matching to the theory with 6 flavors, 
we get 

/ x foi s (M z ) - 0.1184 \ f M t \ 

g 8 (M t ) = 1.1645 + 0.0031 f sl J - 0.00046 - 173.15 1 . (59) 

We determine the MS top-quark Yukawa coupling (y t ) starting from the top-quark pole 
mass (M*) determined from experiments. Averaging measurements from Tevatron and LHC 
experiments, 

( 173.2 ± 0.9 GeV Tevatron [28] 

M t = I 172.6 ± 0.6 ± 1.2 GeV CMS fxj [29] (60) 
{ 174.5 ± 0.6 ± 2.3 GeV ATLAS £j [30], 

we get 

M t = (173.1 ± 0.7) GeV. (61) 

In order to translate this value into a determination of yt we apply: 1) QCD threshold cor- 
rections up to O(c^) [31,32]; 2) complete one-loop electroweak corrections from ref. [33]; 3) 
two-loop 0(aa s ) corrections from ref. [34], including the 0(aa s ) terms due to the renor- 
malization of the Fermi coupling (see sect. 2). As a result, we find, for the MS top Yukawa 
coupling renormalized at the top pole mass M t : 

, , ^ f M t \ ( M h 

y t (M t ) = 0.93587 ± 0.00557 - 173.15 J - 0.00003 f-^^ - 125 

fa s (M z ) - 0.1184 \ . . 

0.00041 f ^ ooQ7 J ± 0.00200 th . (62) 
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Figure 4: The instability scale Aj at which the SM potential becomes negative as a function of the 
Higgs mass (left) and of the top mass (right). The theoretical error is not shown and corresponds 
to a ±1 GeV uncertainty in Mh- 

The 0(aa s ) term, that is the parametrically smallest correction, is equivalent to a tiny shift 
in M t below 0.1 GeV. This effect is well below the O(Aqcd) irreducible non-perturbative 
uncertainty on the top-quark mass determined at hadron colliders (see e.g. ref. [35]), that 
is responsible for the theoretical error in eq. (62). More explicitly, we estimate an irre- 
ducible theoretical error of ±Aqcd ~ ±0.3 GeV in M t from non-perturbative effects, and an 
additional uncertainty of ±0.15 GeV from missing 0(oq) threshold corrections. 

Next, applying the threshold corrections discussed in section 2, we determine the following 
value for the Higgs self coupling in the MS scheme renormalized at the pole top mass: 

A(M t ) = 0.12577 + 0.00205 ~ 125 J - 0.00004 - 173.15 J ± 0.00140 th . (63) 

The residual theoretical uncertainty, that is equivalent to an error of ±0.7 GeV in Mh, has 
been estimated varying the low-energy matching scale for A between Mz and 2M t . 

For completeness, we also include in the one- and two-loop RG equation the contributions 
of the small bottom and tau Yukawa couplings, as computed from the MS 6-quark mass, 
m b (m b ) = 4.2 GeV, and from M T = 1.777 GeV. 
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Figure 5: Regions of absolute stability, meta- stability and instability of the SM vacuum in the Mf- 
M/j plane. Right; Zoom in the region of the preferred experimental range of Mh and Mt (the 
gray areas denote the allowed region at 1, 2, and 3a). The three boundaries lines correspond to 
a s (Mz) = 0.1184 db 0.0007, and the grading of the colors indicates the size of the theoretical error. 
The dotted contour-lines show the instability scale A in GeV assuming a s (Mz) = 0.1184. 

3.3 Phase diagram of the SM 

The final result for the condition of absolute stability is presented in eq. (2). The central 
value of the stability bound at NNLO on Mh is shifted with respect to NLO computations 
(where the matching scale is fixed at fi = M t ) by about +0.5 GeV, whose main contributions 
can be decomposed as follows: 

+ 0.6 GeV due to the QCD threshold corrections to A (in agreement with [14]); 
+ 0.2 GeV due to the Yukawa threshold corrections to A; 

— 0.2 GeV from RG equation at 3 loops (from [12, 13]); 

— 0.1 GeV from the effective potential at 2 loops. 

As a result of these corrections, the instability scale is lowered by a factor ~ 2, for Mh ~ 125 
GeV, after including NNLO effects. The value of the instability scale is shown in fig. 4. 

The phase diagram of the SM Higgs potential is shown in fig. 5 in the M t -Mh plane, 
taking into account the values for Mh favored by ATLAS and CMS data [1,2]. The left 
plot illustrates the remarkable coincidence for which the SM appears to live right at the 
border between the stability and instability regions. As can be inferred from the right plot, 
which zooms into the relevant region, there is significant preference for meta-stability of the 
SM potential. By taking into account all uncertainties, we find that the stability region is 
disfavored by present data by 2a. For M h < 126 GeV, stability up to the Planck mass is 
excluded at 98% C.L. (one sided). 
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Type of error 


Estimate of the error 


Impact on 


M t 


experimental uncertainty in M t 


±1.4 GeV 


a s 


experimental uncertainty in a s 


±0.5 GeV 


Experiment 


Total combined in quadrature 


±1.5 GeV 


A 


scale variation in A 


±0.7 GeV 


Ut 


O(Aqcd) correction to M t 


±0.6 GeV 


Vt 


QCD threshold at 4 loops 


±0.3 GeV 


RGE 


EW at 3 loops + QCD at 4 loops 


±0.2 GeV 


Theory 


Total combined in quadrature 


±1.0 GeV 



Table 1: Dominant sources of experimental and theoretical errors in the computation of the SM 
stability bound on the Higgs mass, eq. (2). 

The dominant uncertainties in the evaluation of the minimum Mh value ensuring absolute 
vacuum stability within the SM are summarized in Table 1. The dominant uncertainty is 
experimental and comes mostly from the measurement of M t . Although experiments at 
the LHC are expected to improve the determination of M t , the error on the top mass will 
remain as the largest source of uncertainty. If no new physics other than the Higgs boson 
is discovered at the LHC, the peculiarity of having found that the SM parameters lie at the 
critical border between stability and metastability regions provides a valid motivation for 
improved top quark mass measurements, possibly at a linear collider. 

The dominant theoretical uncertainty, while reduced by about a factor of 3 with the 
present work, is still related to threshold corrections to the Higgs coupling A at the weak 
scale. Another sizable theoretical uncertainty comes from the fact that the pole top mass 
determined at hadron colliders suffers from O(Aqcd) non-perturbative uncertainties [35]. A 
possibility to overcome this problem and, at the same time, to improve the experimental 
error on M t , would be a direct determination of the MS top-quark running mass from ex- 
periments, for instance from the ti cross-section at a future e + e~ collider operating above 
the ti threshold. In this respect, such a collider could become crucial for establishing the 
structure of the vacuum and the ultimate fate of our universe. 

As far as the RG equations are concerned, the error of ±0.2 GeV is a conservative 
estimate, based on the parametric size of the missing terms. The smallness of this error, 
compared to the uncertainty due to threshold corrections, can be understood by the smallness 
of all the couplings at high scales: four-loop terms in the RG equations do not compete with 
finite tree-loop corrections close to the electroweak scale, where the strong and the top-quark 
Yukawa coupling are large. 

The LHC will be able to measure the Higgs mass with an accuracy of about 100-200 
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Figure 6: Left: The Higgs quartic coupling A at the Planck scale as a function of Mh, with ±3cr 
variations in Mt and a s . Right: Curves in the Mf-M^ plane corresponding to the conditions 
A(Mpi) = (red), f3\(Mp\) = (blue), and to the Veltman's condition of vanishing supertrace 
(green, see text). In all cases the dashed and dotted lines denote the ±3a variation in a s . The gray 
areas denote the experimentally allowed region at 1, 2, and 3a. 

MeV, which is far better than the theoretical error with which we are able to determine the 
condition of absolute stability. 



4 Implications 

4.1 Boundary conditions at the Planck scale 

It is certainly a remarkable coincidence that both A and its beta function /3a nearly vanish 
around the Planck scale. This motivates us to explore in more detail the boundary conditions 
at Mpi required to reproduce the measured values of the SM parameters. In fig. 6 (left) we 
show the prediction for the Higgs quartic coupling A at Mpi, with 3cr bands describing the 
errors in M t and a s . For Mh in the range preferred by LHC, A(Mpi) = can be obtained 
only if M t is ~ 2 GeV below its present central value (i.e. M t < 171 GeV). However, it 
should be noted that A = is neither a fixed point nor a point of enhanced symmetry and 
thus satisfying exactly this condition is not especially meaningful. The best-fit value for A 
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at the Planck scale is small and negative, 
A(M H ) = -0.0144 + 0.0028 ( — - 




- 125 ± 0.0047 Mt ± 0.0018 Qs ± 0.0028 th 



(64) 



where the uncertainties refer to the present la errors in M t and a s , and to the theoretical 
error. 

Equation (64) could be the result of matching the SM at Mpi with a theory in which A 
vanishes at tree level and receives one-loop threshold corrections. However, the contribution 
to high-energy thresholds from the SM couplings at Mpi is typically a few times 10~ 3 . The 
required effect of size 10~ 2 can be obtained from loops of SM couplings only if these involve 
particles with large multiplicities, or else from loop with new (large) coupling constants. As 
evident from fig. 1, the boundary condition of A at high energy is fairly independent of the 
precise value at which we impose it. 

The right plot in fig. 6 shows the dependence on Mh and M t of the boundary conditions 
A(M P1 ) = (red line) and (3\(Mp\) = (blue line). While A = weakly depends on the 
scale at which is evaluated, a more pronounced dependence affects the condition (3\ = (see 
fig. 3). This is because (3\ depends not only on A, but also on other couplings (top Yukawa 
and gauge) that run in the high-energy region. As a result, although (3\(M-p\) = cannot be 
exactly satisfied, the beta function vanishes at scales very near the Planck mass. For instance, 
for M t = 171.0 GeV and Mh = 125 GeV, both A and (3\ simultaneously vanish when they 
are evaluated at a scale equal to 3 x 10 17 GeV. In ref. [16] it was argued that A(Mpi) ps 
and (3\(Mp{) ps could be justified in the case of an asymptotically safe gravitational theory. 

Just for illustration, in the right plot in fig. 6 we also show the Higgs mass implied 
by Veltman's condition [36] that the supertrace of the squared masses of all SM particles 
vanishes at a given scale, here chosen to be Mp\. STrA4 2 (Mpi) = 0. We remark, however, that 
this condition does not carry special information about the power divergences of the theory, 
which are dominated by UV effects, while the supertrace includes only the contribution from 
the IR degrees of freedom in the SM. At any rate, the possibility of a very special fine-tuning 
involving only the SM loop contributions and leading to STr.M 2 (M P i) = 0, implies a Higgs 
mass Mh ~ (135 ±2.5) GeV, which is excluded at more than 3a. Lowering the scale at which 
the supertrace condition is evaluated makes the disagreement even stronger. 

4.2 Higgs inflation from non-minimal coupling to gravity 

The extrapolation of the SM up to very high energy has led to some speculations about the 
possibility of interpreting the Higgs boson as the inflaton. One scenario for Higgs inflation [37] 
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exploits a large non-minimal coupling between the Higgs bilinear and the Ricci scalar R, 
with an interaction Lagrangian £\H\ 2 R. The effect of this interaction is to flatten the Higgs 
potential (or any other potential) above the scale Mpi/ a/£, providing a platform for slow-roll 
inflation. A correct normalization of the spectrum of primordial fluctuations fixes the value 
of the coupling constant £. Using the tree-level potential, one finds £ ~ 5 x 10 4 VX 

This inflationary scenario, attractive for its minimality, suffers from a serious drawback. 
Perturbative unitarity is violated at the scale Mpi/£, signaling the presence of new physical 
phenomena associated with strong dynamics. It is naturally expected that these phenomena 
will affect the scalar potential above Mpi/£ in an uncontrollable way [38]. One solution is to 
add new degrees of freedom that restore perturbative unitarity [39], although the minimality 
of the model is then lost. The procedure advocated by the proponents of this scenario is to 
assume that the strong dynamics will preserve intact the shape of the SM potential, even 
above Mpi/£. Although we find this assumption questionable, it is still interesting to address 
the issue of whether the Higgs data are compatible with this scenario. 

A two- loop analysis of Higgs ^-inflation was developed in ref. [40]. The ordinary SM 
evolution is perfectly adequate below the scale Mpi/£, while the new interaction can affect the 
scalar potential at very high energy. The renormalization procedure above the inflationary 
scale is not unambiguous as, for instance, the renormalization scale differs in the Jordan and 
Einstein frames. Luckily, the slow running of A at high energy makes these issues irrelevant, 
from a practical point of view. A simple SM calculation of the potential is perfectly adequate 
to describe the situation of Higgs ^-inflation (see also the discussion in ref. [14]). 

In practice, the result is that Higgs ^-inflation requires stability of the potential up to 
the inflationary scale Mpi/ a/£- As we are interested in the minimum value of the Higgs 
mass that satisfies this condition, the coupling A at the relevant scale is very small and thus 
the coupling £ is not particularly large, £ < O(10 3 ). Therefore, the resulting restriction is 
stability, as given by eq. (2). If the LHC indication for Mh = 125-126 GeV is confirmed, 
the simplest version of Higgs inflation is disfavored, unless the top mass is about 2a below 
its present central value. However, given the proximity of A(Mpi) to the critical value for 
stability, unknown one-loop threshold corrections near the Planck mass could be sufficient 
to rescue the proposal. It is also interesting that the introduction of a single scalar field at 
the scale Mpi/£ could simultaneously restore perturbative unitarity and cure the potential 
instability [41]. 
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4.3 Higgs inflation from false vacuum 



Alternative proposals for Higgs inflation employ the peculiarity of the SM scalar potential to 
develop a second minimum at large Higgs field values for a very special choice of parameters 
[5, 15]. The possibility of using this new minimum for inflation was first contemplated in 
ref. [6], finding that it implies a viable prediction for the Higgs mass, but also a wrong 
prediction for the amplitude of density fluctuations. The latter result can be cured in non- 
minimal inflationary setups [18,19] without affecting the prediction for the Higgs mass, which 
we now precisely compute. 



The first derivative of the Higgs potential V = X e s(h)h 4 /4: is 



dV 
~dh 



(65) 



Here X e g(h) is the effective coupling defined in eq. (52) and /3 e ff = dX e s/d\ia h. If A c g- becomes 



sufficiently small, the potential can develop a minimum at h = /i mm , such that 



A e ff + 



cff 



o . 



(66) 



h=h„ 



This situation can occur in the proximity of a field value h* where (3 e s vanishes. In the 
neighborhood of h*, we can approximate X e s{h) as 

h 



X eS (h) « A* + 6 In 



h. 



(67) 



where A* is the minimum value of A e fj, such that /3 c fr(A*) = 0. The zero of the (3 function 
insures that the leading log is absent in eq. (67) and thus 6 is a typical two- loop coefficient. 
For the relevant values of the SM parameters, we find 6 = 0.4/(47r) 4 . We are interested 
in a situation in which the field configuration corresponds to a local minimum (while the 
EW vacuum remains the global minimum) and thus we want both A* and b to be positive. 
Using the expansion in eq. (67), we can compute h min and the minimum of the potential 

^niin = ^(^min); 



K exp 



V m 



b 



1 



(6* 



(69) 



The minimum /i mm exists only for extremely small values of the Higgs quartic coupling, 
A* < 6/16. As we vary A* within its allowed range (0 < A* < 6/16) we find that h min is 
always near (e -1 / 4 < /i m i n //i* < 1), while V min can change widely (0 < V min < 6/i min /32). 
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Figure 7: Two-loop SM Higgs potential for Mh = 124, 125, 126 GeV around the critical top mass 
that gives a second minimum around the Planck scale. The various curves correspond to variations 
in M t by 0.1 MeV. 

If the Higgs field is trapped in the false vacuum during the early universe, it can cause in- 
flation. The normalization of the spectrum of primordial perturbations, which is determined 
by Vmin, can be appropriately selected by tuning the ratio A*/6. The main difficulty of this 
scenario is to achieve a graceful exit from the inflationary phase. Two mechanisms have been 
proposed. The first one [18] employs a new scalar field, non-minimally coupled to gravity, 
that slows down the expansion rate, thus allowing for quantum tunneling of the Higgs out 
of the false vacuum. The second mechanism [19] uses a scalar field weakly coupled to the 
Higgs which, during the cosmological evolution, removes the barrier in the Higgs potential in 
a process analogous to hybrid inflation. So, in practice, the minimality of the SM is lost and 
one may wonder if there is any conceptual gain with respect to adding a new scalar playing 
the role of the inflaton. Nevertheless, it is interesting to investigate whether the Higgs and 
top masses are compatible with the intriguing possibility of a false vacuum at large field 
value. 

In practice the above equations amount to saying that the conditions for the existence of 
a second (unstable) vacuum are that A e g- (essentially) vanishes at the same scale at which fl e g 
vanishes. This corresponds to the intersection between the red band (condition A ~ A c fr = 0) 
and the blue band (condition f3\ = 0) in fig. 6 (right). It is remarkable that the SM can achieve 
these conditions, although they require a top mass about 2a below the central value. The 
resulting relation between M h and M t corresponds to the equality in eq. (2), and is precisely 
studied in fig. 7 where we compute for Mh = {124, 125, 126} GeV the predicted top mass and 
show the shape of the potential around the false vacuum. The value of V m i n can be changed 
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Figure 8: NNLO prediction for the Higgs mass Mh in High-Scale Supersymmetry (blue, lower) and 
Split Supersymmetry (red, upper) for tan j3 = {1,2,4,50}. The thickness of the lower boundary at 
tan (3 = 1 and of the upper boundary at tan (3 = 50 shows the uncertainty due to the present la 
error on a s (black band) and on the top mass (larger colored band). 

by tuning A* or, in other words, by accurate variations of Mh and M t . The existence of 
the false vacuum depends critically on the exact values of the SM parameters and requires 
dialing Mh and M t by one part in 10 6 . However, the exact value of the needed top mass has a 
theoretical uncertainty, reduced down to ±0.5 GeV thanks to our higher-order computation. 
Note from fig. 7 that the field value where the false vacuum is positioned is larger than what 
was reported in [6,18]. The corrections in eq. (52) [3,5] are mostly responsible for the larger 
field values found in our analysis. 

4.4 Supersymmetry 

Our higher order computation of the relation between the Higgs mass and the Higgs quartic 
coupling A has implications for any model that can predict A. If supersymmetry is present 
at some scale m, then in the minimal model one finds the tree-level relation 




(70) 
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A dedicated analysis of the resulting prediction for the Higgs mass as function of rh and 
of tan/3 was performed in [42] (see also [43]). We here update the results, including the 
new correction which increases the predicted Higgs mass by an amount that changes with 
M/j. Once more, the main impact of our calculation is the reduction of the theoretical 
uncertainty from ±3 GeV down to ±1 GeV. As a consequence, supersymmetry broken at 
the Planck scale, which requires A > and thus the stability condition eq. (2), is disfavored at 
2a, unless thresholds at Mp\ (or non-minimal couplings) account for the small negative value 
of A(Mpi), see eq. (64). Thresholds at one loop have been computed in [42] including finite 
parts and the constant factor due to the fact that the MS renormalization scheme breaks 
supersymmetry. A simplified, but illustrative, formula for the supersymmetric threshold 
corrections is obtained by taking a common mass M 2 for weak gauginos and higgsinos, and 
a common mass m 4 ~ for the stops, 



SX(M m ) ^-0.25 + 0.121n^ +0.05 In , (71) 

where uia is the mass of the Higgs pseudoscalar. The absence of a scale dependence in 
eq. (71) is a consequence of the approximate cancellation of (3\ around Mpi. Negative values 
of the boundary condition for A thus require stops lighter than higgsinos, winos, and bino. By 
using the full formula and allowing for each sparticle mass to vary by one order of magnitude 
above or below the average mass rh we find 

- 0.006 < A(Mpi) < 0.002 . (72) 

This is insufficient to reach the central value of eq. (64) and thus indicates that supersymme- 
try has to be broken at a scale below Mp\, if the Higgs mass has to match a supersymmetric 
boundary condition. 

Our predictions for the Higgs mass as a function of the supersymmetry breaking scale 
rh are illustrated in fig. 8, in the case of High-Scale Supersymmetry (all supersymmetric 
particles with masses equal to rh) and Split Supersymmetry [44] (supersymmetric scalars 
with masses equal to rh and supersymmetric fermions at the weak scale). We refer the 
reader to ref. [42] for details about the underlying assumptions of the calculation. 

Figure 8 shows not only how Mh ~ 125 GeV disfavors supersymmetry broken at a very 
high scale, but also the well know fact that the usual scenario of weak-scale supersymmetry 
can account for the Higgs mass only for extreme values of the parameters (such as large 
tan/3, heavy stops, maximal stop mixing). In the case of Split Supersymmetry, large values 
of rh are clearly excluded by LHC Higgs searches. On the other hand, (mildly "unnatural") 
scenarios in which the masses of supersymmetric scalars are one-loop larger than the weak 
scale [45] are in good agreement with data. 
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5 Conclusions 



One of the most important questions addressed by the LHC is naturalness. Not only will 
the answer affect our understanding of the mechanism for EW breaking, but it will also 
determine our strategy for future directions in theoretical physics. On one side we have 
the avenue of larger symmetries unifying physical laws in a single fundamental principle; 
on the other side we have new kinds of paradigms, where parameters are not understood 
by naturalness arguments in the context of well-defined effective theories. At present, the 
multiverse is one of the most intriguing options to pursue the latter path. 

If the LHC finds Higgs couplings deviating from the SM prediction and new degrees of 
freedom at the TeV scale, then the most important question will be to see if a consistent 
and natural (in the technical sense) explanation of EW breaking emerges from experimental 
data. But if the LHC discovers that the Higgs boson is not accompanied by any new physics, 
then it will be much harder for theorists to unveil the underlying organizing principles of 
nature. The multiverse, although being a stimulating physical concept, is discouragingly 
difficult to test from an empirical point of view. The measurement of the Higgs mass may 
provide a precious handle to gather some indirect information. 

Once we extrapolate the SM to very short distances, we find that the values of the 
Higgs mass, hinted by the first LHC results (125-126 GeV), lie right at the edge between 
EW stability and instability regions, see fig. 5. Moreover, the slow running of the quartic 
Higgs coupling A in the high-energy regime implies that the instability scale most critically 
depends on the Higgs mass Mh- A small change in Mh (and M t ) can drastically modify our 
conclusions regarding vacuum stability. This special situation motivated us to perform a 
NNLO calculation of the Higgs potential in the SM, which is the subject of this paper. 

Our calculation includes three-loop running for gauge, top Yukawa, and Higgs quartic 
coupling and two-loop matching conditions keeping the leading effects in a s , yt and A. In 
particular, we have computed for the first time the two-loop threshold correction to A, which 
was the most sizable missing ingredient of the NNLO result. The completion of the NNLO 
calculation allows us to reduce the theoretical error in the stability limit on the Higgs mass 
from 3 GeV to 1 GeV. Our final result is shown in eq. (2). After our calculation, the largest 
source of uncertainty comes from the parametric dependence on the top quark mass, which 
leads to a la error of 1.4 GeV in the critical Higgs mass. Of course our calculation reliably 
accounts for IR effects, but ignores possible new (unknown) UV threshold effects occurring 
near the Planck scale. Since our goal is to learn information about physics at very short 
distances, the high-energy corrections constitute an opportunity, rather than a deficiency in 
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the calculation. 

The first lesson that we learn from the SM extrapolation is that the Higgs mass hinted by 
LHC results corresponds to A ~ and /3\ ~ at high energies. This, by itself, is an intriguing 
result because A = is the critical value for stability and it may hide some information about 
Planckian physics. With our precise calculation, we can investigate further the situation. 
We find that, for M h = 125 GeV, \{M Pl ) = -0.014±0.006, see eq. (64). The exact vanishing 
of A(Mpi) is excluded at 2a. Moreover, the smallness of 0\ at high energy implies that A 
remains negative in a relatively large energy range. Indeed, we find that, for Mh = 125 GeV, 
the instability scale develops at 10 11=tl GeV. Quantum tunneling is sufficiently slow to ensure 
at least metastability of the EW vacuum, see fig. 5. The small value of A(Mpi) may indicate 
a radiative origin, although typical one-loop effects of SM couplings appears to be insufficient 
to account for it. We have also shown that /3\ varies more rapidly at high energy and vanishes 
at a scale of about 3 x 10 17 GeV. 

The stability of the SM potential is a crucial issue for models of inflation that employ 
the Higgs boson. We have analyzed several proposals showing that present data disfavor 
them at 98% C.L. These models can still be viable if the top quark mass turns out to be less 
than about 171 GeV or if new physics around Mp\ slightly modifies the shape of the Higgs 
potential. The latter possibility, although fairly plausible, limits the predictability and the 
minimality of the approach. We have also updated previous predictions for the Higgs mass 
in High-Scale Supersymmetry and Split Super symmetry. 

It is natural to try to speculate on the possible meaning of the near vanishing of A and 
(3\ around the Planck scale. The coupling A = is the critical value that separates the 
ordinary EW phase from a phase in which the Higgs field slides to very large values. It is 
noteworthy that the hierarchy problem can also be interpreted as a sign of near criticality 
between two phases [46]. The coefficient m 2 of the Higgs bilinear in the scalar potential is 
the order parameter that describes the transition between the symmetric phase (m 2 > 0) 
and the broken phase (m 2 < 0). In principle, m 2 could take any value between — Mpj and 
+Mpj, but quantum corrections push m 2 away from zero towards one of the two end points 
of the allowed range. The hierarchy problem is the observation that in our universe the 
value of m 2 is approximately zero or, in other words, sits near the boundary between the 
symmetric and broken phases. Therefore, if the LHC result is confirmed, we must conclude 
that both m 2 and A, the two parameters of the Higgs potential, happen to be near critical 
lines that separate the EW phase from a different (and inhospitable) phase of the SM. We 
do not know if this peculiar quasi-criticality of the Higgs parameters is just a capricious 
numerical coincidence or the herald of some hidden truth. 
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The occurrence of criticality could be the consequence of symmetry. For instance, su- 
persymmetry implies m 2 = 0. If supersymmetry is marginally broken, m 2 would remain 
near zero, solving the hierarchy problem. But if no new physics is discovered at the LHC, we 
should turn away from symmetry and look elsewhere for an explanation of the near-criticality 
of m 2 . 

The critical value A = could be justified by symmetry reasons. For instance, if the Higgs 
is a Goldstone boson, its potential vanishes and both m 2 and A are zero. The non-vanishing 
top Yukawa coupling prevents this possibility to be realized exactly. Radiative corrections 
then completely spoil the solution to the hierarchy problem, but could generate a small and 
negative value of A at the Planck scale, compatible with our results. Supersymmetry broken 
at high scales could also account for a vanishing boundary condition of A, if tan/3 = 1. 
However, we have shown that, unless the Higgs has strong couplings with new states that 
live at Planckian energy, the supersymmetric boundary condition cannot be satisfied at Mpi, 
see fig. 8. 

Alternatively, criticality could be the consequence of dynamics. If transplanckian dynam- 
ics induce a large anomalous dimension for the Higgs field, the matching condition at Mpi 
of the quartic coupling A could be very small, while the top Yukawa coupling remains siz- 
able. It was argued in ref. [16] that gravity itself could be responsible for a large anomalous 
dimension of the Higgs in the transplanckian region. 

It is known that statistical systems often approach critical behaviors as a consequence 
of some internal dynamics or are attracted to the critical point by the phenomenon of self- 
organized criticality [47]. As long as no new physics is discovered, the lack of evidence for 
a symmetry explanation of the hierarchy problem will stimulate the search for alternative 
solutions. The observation that both parameters in the Higgs potential are quasi-critical 
may be viewed as evidence for an underlying statistical system that approaches criticality. 
The multiverse is the most natural candidate to play the role of the underlying statistical 
system for SM parameters. If this vision is correct, it will lead to a new interpretation of our 
status in the multiverse: our universe is not a special element of the multiverse where the 
parameters have the peculiarity of allowing for life, but rather our universe is one of the most 
common products of the multiverse because it lies near an attractor critical point. In other 
words, the parameter distribution in the multiverse, instead of being flat or described by 
simple power laws (as usually assumed) could be highly peaked around critical lines because 
of some internal dynamics. Rather than being selected by anthropic reasons, our universe is 
simply a very generic specimen in the multitude of the multiverse. 
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The indication for a Higgs mass in the range 125-126 GeV is the most important result 
from the LHC so far. If no new physics at the TeV scale is discovered, it will remain as 
one of the few and precious handles for us to understand the governing principles of nature. 
The apparent near criticality of the Higgs parameters may then contain information about 
physics at the deepest level. 
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Appendix: SM Effective Potential up to two-loops 

The SM effective potential is known up to two-loops [20]. We present here its explicit 
expression in the MS scheme and the Landau gauge. The tree-level part is (in this section 
we denote the Higgs field by <fi to avoid confusion with the symbol h defined later) 

K,(0) = -im 2 2 + |A0 4 . (73) 

The one-loop Coleman- Weinberg potential [48] is 

V x {(j>) = | [-12m$(L t - 3/2) + 6m 4 w {L w - 5/6) + 3m|(L 2 - 5/6) 

+mi(L h - 3/2) + 3m 4 x (L x - 3/2)] (74) 

where k = l/(167r 2 ), m 2 = y 2 (j) 2 /2 is the top mass squared, m\ = —m 2 + 3\<j) 2 the Higgs 
mass squared, rn 2 {4>) = —m 2 + \<p 2 is the Goldstone mass squared and L t = \n(m 2 /Q 2 ), etc. 

We split the two-loop potential in different pieces according to their diagrammatic origin. 
We use the short-hand notation t = m 2 , w = m 2 w , z = m 2 ., h = ml, \ = m x and we neglect 
the bottom Yukawa coupling. The important top Yukawa contribution is 

V Y = ^y 2 t K 2 {2J tt -AJ tg -2J th + {At-h)I tth + 2{t-x)h g v-xhtz] ■ (75) 
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There is a purely scalar piece 



a purely gauge part 



3 

= ~K 2 X [hJgg + 2J hg + J hh ~ A\4> 2 (l hgg + I hhh )] 



(76) 
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a fermion-gauge boson part 3 (which includes the important QCD piece) 
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in ref. [20] contains an erratum for this piece, with an extra factor 3 for the 
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3 Noticc that the first paper in 
lepton-lepton-^ contributions. 
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and a scalar-gauge boson part 
"1 
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The functions A, J and / are 

A x = A[x] = x(L x - 1) , 
J xy = J[x,y] = A[x]A[y] , 

I xyz = I[x,y,z] = -[(x-y-z)L y L z + (-x + y-z)L x L z + (-x-y + z)L x L y ] 
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where L x = \n(x/Q 2 ) and 
£[x,y,z] = R 
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where i? 2 = x 2 + y 2 + z 2 — 2xy — 2xz — 2y^ and Li 2 (x) is the dilogarithm function. The 
above expression is valid for R 2 > 0, while for R 2 < the analytical continuation should be 
used instead. Some particular cases of the previous functions which are useful to evaluate 
the potential are 
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The 2-loop subleading corrections to A e fr not included in eq. (52), obtained from the full 
2-loop correction (but setting A — > 0, as this coupling is particularly small at high field 
values) are explicitly given by 
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where f 11{BW = l,x/y), 

r p = lnfKpe 27 ^] , r t / TO = lnf/^//^] , r( t -w)/t = bi[(«t - ««,)/««,] , 

and so on. 
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